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ABSTRACT 



Context. A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared 
to the predictions in the framework of the "standard" cosmological model. This is known as the "arc statistics problem". Recently, few other 
inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the 
^> high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may 
■^j" ■ be gravitational lenses stronger than expected. 

' Aims. We aim at better understanding these issues by comparing the lensing properties of well defined cluster samples with those of a large set 
C^> of numerically simulated objects. 

■ Methods. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray 
| flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample 
f^*) , of 12 MACS clusters at z > 0.5 presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). 

t-H ' Results. We find that three clusters in the MACS sample have lensing cross sections and Einstein ring sizes larger than any simulated cluster 
t-H i in the MareNostrum Universe. We use the lensing cross sections of both simulated and real clusters to estimate the number of giant arcs that 
^) . should arise by lensing sources at z = 2. We find that simulated clusters produce ~ 50% less arcs than observed clusters do. The medians of 
• i-H ' the distributions of the Einstein ring sizes differ by ~ 25% between simulations and observations. We estimate that, due to cluster triaxiality 
, and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the 
S— i lensing analysis should be up to a factor of ~ 2 larger than expected from the ACDM model. In particular, we predict that for ~ 20% of the 
clusters in the MACS sample the lensing-derived concentrations should be higher than expected by more than ~ 40%. 

Conclusions. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the 
basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, 
the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending 
our comparison. 

Key words. Keywords should be given 



1. Introduction 

Strong lensing is a widely used method fo r invest i gatin 



the inner struc t ure of galaxy clus t ers (s e e iKovnerl I1989I 



Bergmann et al. Il990t Kneib et alj l2003b Broadhurst et al 
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2005 



Cacciato et al.l 2006; Li esenborgs et al. 2009; ICoe et al 



2010, for some examples). Moreover, the statistical analysis of 
strong lensing eve nts in clusters was also proposed as a tool 
for cosmology (e.g. lBartelmann et al.l ll998). The abundance of 
strong lensing events, such as gravitational arcs and multiple 
images, is expected to be higher in cosmological models where 
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the growth of the cosmic structures is faster at earlier epochs, 
such as models where some dynamical dark energy starts 
dominating the expansion of the universe at earlier epoch s 
dBartelmann et al.l2003l;lMacci6l200"5l:lMeneghetti et alj2005l) . 
Thus, in these models a larger number of potential lenses pop- 
ulates the universe up to high redshift (z > 0.5). The clus- 
ter concentrations are found in numerical simulations to re- 
flect the density of the universe at the cluster epoch of for- 
mation. Clusters f orming earlier have higher concentrations 
dDolag et al.l2004l) and are expected to be more efficient lenses. 

In the era precision cosmology, strong lensing statistics 
cannot be as competitive as other cosmological probes for con- 
straining cosmological parameters. These are all supporting 
the "concordance model" ACDM mode l, which became the 
standard scenario of structure formation (Komatsu et al. 2009; 
Ries s et all 19981 12004 IPerlmutter et al.lll999l:lEisenstein et al.1 
2005; Percival et al.1 120071) . However, previous attempts of 
using strong lensing statistics as a cosmo logical tool have 



produc ed controversial results. In particular, Bartelma nn et al 



( 1998b . studying the lensing properties of a set of numerically 
simulated galaxy clusters, argued that the ACDM cosmologi- 
cal model fails at reproducing the observed abundance of giant 
gravitational arcs by almost an order of magnitude. This incon- 
sistency between strong lensing and other observational data is 
known as the "arc statistics problem". Its nature is still under 
debate. A long seri es of papers have tried to falsify the theoret- 
ical predictions of Bartelma nn et al.l dl998l) . Although several 
limits where found in their simulations, which could not prop- 
erly capture sever al important features of both the lenses and 
the sources (e.g. iDalalet all |2004 IWambsgansset all 1 20041; 
Meneghetti et alJ2003l lTorri et al.ll2004 IPuchwein et aLfcOOa 
Meneghe tti et al.1 120071; IWambsganss et alj 120081; iMead etaD 
201 Oh . the controversy is not yet solved. Moreover, it has been 
recently enforced by several other observations of strong lens- 
ing clusters, which seem to indicate that 1) some galaxy clus- 
ters have very extended Einstein rings (i.e. critical lines) whose 
abundances can hardly be reproduced by cluster models in 
the framework of a ACDM c osmology ( Broadhurst & Bar kana 
2008| iTasitsiomi et al.1 12004), and 2) few clusters, for which 



high quality strong and weak lensing data became available, 
have concentratio ns which are way too large compared to 



the expectations (IBroadhurst et al.1 2008; Zi trin et al.l [2009b) 



These evidences push in the same direction of the "arc statis- 
tics" problem, in the sense that they both suggest that observed 
galaxy clusters are too strong lenses compared to numerically 
simulated clusters. 

Understanding the origin of these mismatches between the- 
oretical predictions and observations is fundamental, as these 
may evidence a lack of understanding of the cluster physics, 
which may be not well implemented in the simulations, or, 
conversely, highlight some inconsistencies between the ACDM 
scenario and the properties of the universe on small scales. So 
far a comparison between theoretical predictions and observa- 
tions has been complicated by the lack of systematic arc sur- 
veys but also by the fact that different approaches were used to 
analyze simulations and observations (Mene ghetti et all2008l) . 
In this paper, we propose a novel approach, whose principal 
aim is to eliminate most of the assumptions used in the pre- 



vious works. It consists of analyzing observed galaxy clusters, 
for which detailed mass models are available through strong, 
and possibly also weak lensing observations, fully consistently 
to numerical simulations. The deflection angle maps provided 
by the lens models are used for performing ray-tracing sim- 
ulations and for lensing the same source population used in 
numerical simulations. We attempt a comparison between the 
properties of simulated clusters in the MareNostrum Universe 
and those of a complete sample of X-ray luminous MACS clus- 
ters, for wh i ch stron g lensing models were recently derived by 
Zitri n et al. (l2010ah . 

The plan of the paper is as follows. In Sects. 2 and 3 we 
describe the numerical and the observed cluster samples used 
in this study. In Sect. 4 we illustrate the analysis performed to 
measure the lensing cross sections and the Einstein ring sizes. 
In Sect. 5 we discuss the comparison between simulations and 
observations, describing how we construct mock cluster cat- 
alogs to simulate a MACS-like survey, and showing the sta- 
tistical distributions of both cross sections and Einstein rings. 
Finally, we use the simulated clusters for estimating the con- 
centration bias expected for the MACS sample. Sect. 6 is dedi- 
cated to the discussion and the conclusions. 

2. The MareNostrum Universe 

Our theoretical expectations are based on the analysis of the 
simulated clusters contained in the MareNostrum Universe. A 
detailed descript ion of the analysis perform ed on t hese objects 



can b e found in Meneghetti et al. (2010a) and in Fedeli et al 



(2010). We briefly summarize the most relevant aspects of the 
analysis here. 

The MareNostrum Universe (Gottlo ber & Yepes 2007 ) is 
a large-scale cosmological non-r adiative SPH s imulation per- 
formed with the Gadget2 code (Springel 2005). It was per- 
formed assuming a ACDM cosmological background with 
WMAP1 normalization, namely Q. m .o = 0.3, £2a,o = 0.7 and 
<X8 = 0.9 with a scale invariant primordial power spectrum. 
The simulation consists of a comoving box size of 500 h~ l Mpc 
containing 1024 3 dark matter particles and 1024 3 gas particles. 
The mass of each dark matter particle equals 8.24 x 10 9 M G h~ l , 
and that of each gas particle, for which only adiabatic physics is 
implemented, is 1.45 x 10 9 M Q h~ l . The baryon density param- 
eter is set to Qb,o = 0.045. The spatial force resolution is set to 
an equivalent Plummer gravitational softening of 15 IT 1 kpc, 
and the SPH smoothing length was set to the 40th neighbor to 
each particle. 

As described in 



Meneghetti et al.l (l2010al) . we extracted 



from the cosmological box all the cluster-sized halos, which 
were subsequently analyzed using ray-tracing techniques. Each 
halo was used to produce three different lens planes, obtained 
by projecting the cluster mass distribution along three orthog- 
onal lines of sight. These planes of matter were used to lens 
a population of elliptical sources of fixed equivalent radius of 
0.5" placed on a source plane at redshift z s = 2. This analysis 
was performed on all the halos found between z\ = and z\ — 2. 
The strong lensing clusters were classified in three main cate- 
gories, namely 1) clusters with resolved critical lines, i.e. which 
can produce strong lensing features such as multiple images of 
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the same background source; 2) clusters with non-zero cross 
section for giant arcs, i.e. clusters which are potentially able 
to distort the images of background galaxies such to form arcs 
with length-to-width (L/W) ratios larger than 7.5; 3) what we 
called "super-lenses", i.e. clusters whose lensing cross section 
for giant arcs is larger than 10~ 3 /z~ 2 Mpc 2 . This definition is 
based on simulations in cluding observational noises per formed 
with the SkyLens code dMeneghetti et al.ll2008ll2010bh . which 
show that, observing a cluster of galaxies with such cross sec- 
tion using ~ 3 HST orbits in the /-band, the expected num- 
ber of giant arcs in the cluster field is ~ 1. Taking advan- 
tage of the large size of the cluster sample, we could char- 
acterize statistically the strong lensing cluster population in 
the MareNostrum Universe, correlating the lensing strength 
to several cluster properties, such as their mass, shape, orienta- 
tion, concentration, dynamical state, and X-ray emission. With 
~ 50, 000 strong lensing clusters found in the MareNostrum 
Universe, this sample is the largest ever used for strong lensing 
studies. 

3. Strong-lensing analysis of The MACS 
high-redshift cluster sample 

In a recent paper, Zitrin et al. (2010) showed the results of 
the strong lensing analysis of a sample of 12 very lumi- 
nous X-ray galaxy clusters at z > 0.5 using HST/ACS im- 
ages. This is a complete sample of clusters with X-ray flux 
fx > 1 x 10~ 12 erg s _1 crrT 2 in the 0.1-2.4 keV band, which 
was defined by Ebeling et al. (2007). Later, it was targeted 
in several follow-up studies, including deep X-ray, SZ and 
HST imaging. For example, the detection of a large-scale fila- 
ment has been repor ted in the case of MACS J0717.5+3745 by 
Ebeling et al.l (2004), for which many multip l y-lense d images 



have been recently identified by Zit rin et al.1 (12009al) . reveal- 



ing this object to be the largest known lens with an Einstein 
radius equivalent to 55" (for a source at z ~ 2.5). In the 
case of MACS Jl 149.5+2223 dZitrin & Broadhurstl 120091) . a 
background spiral galaxy at z — 1.49 ( Smith et al. 2009) has 
been shown to be multiply-lensed into several very large im- 
ages, with a total magnification factor of ~ 200. Another 
large multiply-lensed sub-mm source at a redshift of z - 
2.9 has been identified in MACS J0454. 1-0300 (also referred 



to as MS 0451.6-0305; Takataetal. 2003; Borvs et al. 2004 



Berc iano Alba et al.1 120071 \2Q 1 Oh, MACS J0025 .4- 1 222 



found to be a "bullet cluster"-like (Bradacetal. 2008) 



was 
and 

other MACS clusters ha ve been recently u sed for an exten- 
sive arc statistics study ( Horesh et al. 2010h . The X-ray data 
available for this sample (see Ebeling et al. 2007) along with 
the high-resolution HST /ACS imaging and additional SZ data 
(e.g. lLaRoque et alJEooll) make these 12 high-redshift MACS 
targets particularly useful for understanding the nature of the 
most massive clusters. 

The strong-lensing modeling of this sample, published in 
full in Zitrin et al. (2010) and summarized briefly here, is mo- 
tivated by the successful minimalistic approach of Broadhurst 
et al. (2005) to lens modeling, simplified further by Zitrin et al. 
(2009b). This simple modeling method relies on the assump- 
tion that mass traces light so that the galaxy distribution is 



the starting point of the mass model, and additional flexibil- 
ity between the dark matter and galaxies is allowed through the 
implementation of external shears. Still, the method involves 
only 6 free parameters, enabling easier constraints on the mass 
model since the number of constraints has to be equal or larger 
to the number of parameters in order to get a reliable fit. Two 
of these parameters are primarily set to reasonable values so 
only 4 of these parameters have to be constrained initially, 
which sets a very reliable starting-point using obvious systems. 
Recently we have f urther tested this assumption in Abell 1703 
dZitrin etal.ll2010bl), where we also ap plied the non-parametric 
technique of lLiesenborgs et al ] d2006h for comparison. This lat- 
ter technique employs an adaptive grid inversion method and 
does not make any prior assumptions of the mass distribution, 
yielding a very similar mass distribution to our parametric tech- 
nique and hence confirming the assumption that mass gener- 
ally traces light. In addition, it has been found independently 
that SL methods based on parametric modeling are accurate at 
the level of few percents at predicting the projected inner mass 
dMeneghetti et alj|2010bh . 

The mass distribution is therefore primarily well con- 
strained, uncovering many multiple images which can be then 
iteratively incorporated into the model, by using their redshift 
estimation and location in the image-plane (e.g., Abell 1689, 
Broadhurst et al. 2005, C10024, Zitrin et al. 2009a). In the par- 
ticular case of the 12 high-z MACS clusters, in most of the clus- 
ters the multiple-images found or used in Zitrin et al. (2010a) 
currently lack redshift information, so that the mass profiles 
of most of these clusters could not be well constrained but 
only roughly estimated by assuming crude photometric red- 
shifts. This however still allows an accurate determination of 
the critical curves for any given multiply-lensed source, as the 
critical curves (and the mass enclosed within them) are not de- 
pendent on the mass profile and are relatively invariant to the 
model parameters, enabling us to securely compare these prop- 
erties to simulations as done here. In addition, note that we 
compare the Einstein radius for sources at z s - 2. Due to the 
lensing-distance ratios, an over estimation of the source red- 
shift by Az ~ 0.5 would only increase the projected mass and 
the observed Einstein radius for a source at z s = 2, thus result- 
ing only in a growth of the discrepancy between observations 
and the ACDM simulations presented here. On the other hand, 
underestimating a source redshift would in practice decrease 
the observed Einstein radius and projected mass for a source 
at z s =2 by less than 10%, thus insignificantly influencing the 
results. 



4. Analysis 

4.1. Lensing cross sections 

The efficiency of a galaxy cluster in producing arcs with a given 
property can be quantified by means of its lensing cross section. 
This is the area on the source plane where a source must be 
placed in order to be ima ged as an arc with t hat pro perty. 

d2010ah 



As explained in Meneghe tti et al 



and in 

Meneghetti et al.l ([2003), the lensing cross sections are derived 
from the deflection angle maps, which are in turn calculated 
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by m e ans of ray-tracing methods (see also iMeneghetti et al. 
20051: iFedeli et all l2006h . Bundles of light rays are traced 
from the observer position back to the source plane. This is 
populated with an adaptive grid of elliptical sources, whose 
spatial resolution increases towards the caustics, in order to 
artificially increase the number of highly magnified images. 
In the following analysis, a statistical weight, w,-, which is 
related to the spatial resolution of the source grid at the source 
position, is assigned to each source. If a is the area of one 
pixel of the highest resolution source grid, then the area on 
the source plane of which the z'-th source is representative is 
given by A, = aw t . The images are analyzed individually by 
m easuring their lengths an d widt hs using the method outl ined 



Bartelmann et al 



dl998b and in IMeneghetti et al] d2000h . 



We define the lensing cross section for giant arcs, <x, as 



(1) 



where the sum is extended to all the sources producing at least 
one image with L/W > 7.5. 

4.2. Einstein rings 

While the lensing cross section is defined on the source plane, 
the Einstein ring is defined on the lens plane. Although the 
word "ring" is appropriate only in the case of axially symmet- 
ric lenses, several authors have used it to indicate the tangential 
critical line of lenses with arbitrary shapes. This is the line 6 t 
defined by the condition 



Me t )=n;\0 t ) = o, 



(2) 



where A, is the inverse tangential magnification, yT t x . The mag- 
nification is related to the lens convergence, k, and shear, y, via 
the equation 

1 



MO) = 



1-40) -y(0) 



(3) 



For axially symmetric lenses the following relation holds be- 
tween k and y: 



y(0) = k{9) - k{9) 



(4) 



where ~k{9) indicates the mean convergence within a circle of 
radius 9. Using the above formulas, the Einstein ring of an axi- 
ally symmetric lens is defined as the distance 9e from the lens 
center where 



1-k(6e)=0 



(5) 



i.e. as the radius of a circle enclosing a mean convergence of 
1 . We remind that the convergence is related to the lens surface 
density S and to the critical surface density for lensing E cr via 
the equation 

1(0) 



K(0) 



(6) 



Thus, the mean surface density within the Einstein ring of an 
axially symmetric lens is 



AnG D Xs D\ ' 



(7) 



o 
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Fig. 1. The median Einstein rings of a lens with increasingly 
larger ellipticity of the lensing potential. Black, blue, and red 
solid lines show the tangential critical lines of lenses with ellip- 
ticity 0, 0.1, and 0.2, respectively. The corresponding median 
Einstein rings are given by the solid, dotted, and dashed lines 
of the same colors. 



100 




100 



3 E,med 



Fig. 2. The comparison between equivalent and median 
Einstein ring sizes for a sample of clusters with mass M > 
5 x I0 u h~ l M Q at z > 0.5 extracted from the MareNostrum 
Universe. Each cross represents a cluster projection. The 
dashed red line indicate the bisector of the %med - 9 E , eqv plane. 



where the D S ,D\, and Di s denote the angular diameter distances 
between the observer and the source plane, between the ob- 
server and the lens plane, and between the lens and the source 
planes, respectively. 

This is strictly applicable only to axially symmetric lenses, 
but this definition of Einstein ring has b een exported to ar- 
bitrary lenses by s everal authors (see e.g Zitr in et al.l l2010a: 
Richard et alJ l2010). which define the equivalent Einstein ring 
size, % e qv, via Eq. [5] In this paper, we follow a different ap- 
proach. We define a median Einstein ring size, f^med, which is 
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the median distance of the tangential critical points from the 
cluster centetQ : 



%med = median(0,) . 



(8) 



This definition better captures the important effect of shear 
caused by the cluster substructures, whose effect is that of elon- 
gating the tangential critical lines along preferred directions, 
pushing the cri tical points to distances where k is well below 
unity (see also Bartelman n et al. I ll995h . For example, an axi- 
ally symmetric lens embedded in an external shear has a larger 
%med compared to an isolated axially symmetric lens. The 
same argument applies to lenses whose lensing potential is el- 
liptical rather than spherical. This is shown in Fig.Q] where we 
display the tangential critical lines of a lens with NFW density 
profile whose iso-potential contours have ellipticities e - ^r^i 
equal to 0. 0.1, and 0.2. The solid, dotted, and dashed lines 
show the respective median Einstein rings, whose radius is de- 
rived from Eq. [8] 

In Fig. [2] we show a comparison between the two defini- 
tions of Einstein ring size, when they are applied to a sample 
of massive clusters (M > 5 x 10 14 /z _1 M ) at redshift z > 0.5 
taken from the MareNostrum Universe. Each cross represents 
a cluster projection in the % me d - %eqv plane. Clearly, al- 
most all the points lay in the bottom-right part of the dia- 
gram, below the bisector shown by the dashed red line. Thus, 
#E,med ^ #E,eqv for the vast majority of the clusters. This in- 
dicates that these systems are typically non-axially symmet- 
ric an d contain many s ubstructures that enhance their shear 
fields. Meneg hetti et alj d2007l) showed that asymmetries and 
substructures contribute significantly to the lensing cross sec- 
tion for giant arcs. Therefore, we expect that % me d correlates 
much better with cr than 6>E,eqv does. This is shown in Fig. [3] In 
the left and in the right panels we plot the lensing cross section 
vs the equivalent and the median Einstein ring size for the same 
clusters used in Fig. [2] The figure shows that, when the median 
Einstein radius is used, all the data points lay very close to a 
line in the log(cr) - log(0E,med) plane, whose equation is 



log(cr) = (1.79 ± O.O4)log(0 E , med ) - (5.16 ± 0.05) 



(9) 



The Pearson correlation coefficient is r — 0.94, confirming 
that the correlation between the two plotted quantities is very 
strong. 

As shown in the left panel, using the equivalent Einstein 
radius, the scatter is substantially larger and the correlation in 
much worse. The Pearson coefficient is r — 0.75 in this case. 
For this reason, in the following analysis we prefer to use the 
median, rather than the equivalent Einstein radius. 

The tight correlation existing between cr and £>E,med high- 
lights the strong connection between the arc statistics and the 



1 When dealing with mass distributions characterized by asymme- 
tries and substructures, the cluster center may not be easily defined. 
In this work, the cluster center is determined by smoothing the pro- 
jected cluster mass distribution with a gaussian kernel with a FWHM 
of 30 kpc. This is meant to erase the local peaks related to the pres- 
ence of galaxy-scale halos. The center of the cluster is then defined as 
the location of the maximum of the smoothed projected mass distribu- 
tion. This is only one of the possible choices for the cluster center. Of 
course, the same definition holds for simulated and real clusters. 



Einstein ring problems. If the Einstein ring sizes were too large 
for the ACDM model, then we would also observe an excess of 
giant arcs compared to the expectations. 

4.3. Analysis of the MACS sample 

As we pointed out in Sect. Q] our purpose is to perform a con- 
sistent comparison between real and simulated galaxy clusters. 
For this reason, we use the deflection angle maps obtaine d for 
the MACS cluster reconstructions of Zitri n et all (2010a I and 
we use them to perform lensing simulations with the same 
methods used to analyze the clusters in the MareNostrum 
Universe. For each cluster, we measure the lensing cross sec- 
tion and the Einstein radius. Table Q] summarizes the results 
found for the 12 MACS clusters used in this paper. For each 
cluster in the sample, we also list the redshift in the second col- 
umn. Note that all but two clusters have lensing cross sections 
for giant arcs cr > 10~ 3 /z~ 2 Mpc 2 . On the basis of the classifica- 



tion proposed in Meneg hetti et al.l (l2010al) . these clusters could 
be classified as super-lenses. 



Table 1. The lensing cross sections, the median and the equiv- 
alent Einstein radii of the 12 MACS clusters used in this 
work. The values r eported for #E,eqv are taken from Table 2 of 
Zitrin etal.ld2010ah . 



MACS 


z 


cr 


#E,med 


$E,eqv 






[10- 3 r 2 Mpc 2 ] 


[arc sec] 


[arcsec 


J0018.5+1626 


0.545 


3.27 


31.5 


24 


J0025.4-1222 


0.584 


15.27 


55.8 


30 


J0257. 1-2325 


0.505 


7.89 


54.2 


39 


J0454. 1-0300 


0.538 


0.98 


13.9 


13 


J0647.7+7015 


0.591 


3.07 


27.5 


28 


J07 17.5+3745 


0.546 


14.2 


71.1 


55 


J0744.8+3927 


0.698 


2.30 


32.2 


31 


J091 1.2+1746 


0.505 


0.25 


10.5 


11 


Jl 149.5+2223 


0.544 


2.69 


25.0 


27 


J1423.8+2404 


0.543 


1.64 


21.5 


20 


J2129.4-0741 


0.589 


8.48 


46.0 


37 


J2214.9-1359 


0.503 


2.22 


26.0 


23 



It is very interesting to note that the #E,med - cr relation found 
for the MACS clusters is consistent with that measured in the 
MareNostrum Universe. This is shown in Fig. [4] Each dia- 
mond indicates the position of a MACS cluster in the %med - cr 
plane. The red solid line shows the best fit relation given in 
Eq.|9] derived from the clusters in the MareNostrum Universe. 
The MACS clusters nicely follow the same relation found in the 
simulations. 



5. Comparison between simulations and 
observations 

We now proceed at comparing statistically the strong lensing 
properties of the MACS clusters to those of the clusters in the 
MareNostrum Universe. 
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Fig. 3. Lensing cross section for giant arcs vs Einstein ring size for a sample of clusters with mass M > 5 x 10 1 A h~ x M Q at z > 0.5 
extracted from the MareNostrum Universe. The two panels refer to the two different definitions of the Einstein ring radius. On 
the left, we use the equivalent Einstein radius, while on the right we use the median Einstein radius. The red solid line in the right 
panel indicates the best linear fit relation between log(cr) and log(#E,med)- 
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Fig. 4. The <x - #E,med relation for the 12 MACS clusters used 
in this work. Each diamond indicates a cluster, while the 
red line is the best-fit relation found for the clusters in the 
MareNostrum Universe and given in Eq.|9] 



5.1. Correction of the X-ray luminosities 

The MACS sample is an X-ra y flux selected samp le. The se- 
lection method is discussed in lEbeling et al. (2007) and sum- 
marized later in this paper. As explained in Meneghetti et al 



(2010a) we have measured the X-ray emission of all clusters 
in the MareNostrum Universe, which, in principle, would al- 
low us to apply the same selection function used for the MACS 
sample to the simulated clusters. However, the MareNostrum 
Universe is non-radiative simulation and it is well known 
that the X-ray properties of galaxy clusters can be reproduced 
only by means of a more sophisticated description of the gas 
physics. In particular, the X-ray luminosity of low mass clus- 
ters is known to be over-predicted in adiabatic simulations, re- 



sulting in a much shallower M - Lx relation than observed. A 
discussion about how the M - Lx relation changes dependin g 
on the different gas physics can be found in I Short et al. (2010). 
The authors of this paper show that only including cooling and 
some mechanism, like pre-heating or AGN feedback, to heat 
the gas and prevent it from reaching high central densities, 
can the simulation match the slope of the observed M - Lx 
relation at low re dshift, as derived from the REXCESS data 
( Pratt et al1l2009b . They also provide analytic formulas for de- 
scribing the redshift evolution of the X-ray scaling relations in 
different kind of simulations. In particular, the M - Lx relation 
is parametrized as follows: 



E(zT 7/3 L x 



where 



C(z) 



M 



^A,0 



E(z) = ^n m , (l+z) 3 +£^ 
for a spatially flat ACDM cosmological model, and 
C(z) = C (l+zf . 



(10) 



(11) 



(12) 



The parameters Co, a, and /? are estimated by fitting the scal- 
ing relations of the simulated clusters. The best fit values are 
reported in Table 2 and 3 of IShort et al.l (120 1 Ob for different 
gas physics implemented in the simulations. The mass Mq is 
5 x lO 14 ^ 1 M Q . 

Using an high re dshift sample of X-ray clusters from 
Maughan et al.l (12008b . the authors find that the scaling rela- 
tions derived by including an AGN feedback model (FO run) 
evolve broadly consistently with the ob servational data. F rom 
Eq. [10] and the best fit values found by IShort etail d2010h for 
their adiabatic (GO) and FO runs, we derive the correction we 
should apply to the X-ray luminosities of our simulated clus- 
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ters in order to facilitate a comparison to the observations. This 
correction is estimated as 



0.316(1 +z) a991 ( — 
M 



(13) 



5.2. Simulating MACS-like surveys 

We use the MareNostrum Universe as a reference to build 
up mock catalogs of MACS-like clusters. The simulated clus- 
ters are distributed on the sky within a spherical shell between 
Z - 0.5 and z = 0.7. This volume is subdivided into seven sub- 
shells, which are populated with objects taken from the snap- 
shots at z = 0.5, 0.53, 0.56, 0.59, 0.63, 0.66 and 0.68 of the sim- 
ulation. When necessary, we replicate the clusters in the cosmo- 
logical box such to reproduce the expected number of halos in 
each shell. The number of replicates is determined by the ratio 
between the shell volume and the simulation volume. For the 
lensing analysis, each cluster is projected along three different 
lines of sight. Every time we replicate a cluster, we randomly 
choose the line of sight along which it is observed. 

The sample of MACS clusters used in this work was con- 
structed by using the following selection criteria: 

- the X-ray flux in the band [0.1-2.4] keV is f x > 1 x 10~ 12 
erg s _1 cirr 2 ; 

- the clusters are observable from Mauna Kea: \b\ > 20°, 
-40° < 5 < 80°; 

- the redshift is z > 0.5 (the most distant cluster is at z = 
0.698, as shown in Table [TJ- 

Applying the cuts to our mock cluster catalogs, we finally de- 
fine a MACS-like sample for our comparison. We repeat the 
procedure to generate the cluster catalogs 10 times to partially 
account for the cosmic variance. 



5.3. Distributions of lensing cross sections 

A comparison between the distributions of the lensing cross 
sections for giant arcs in the simulated and in the observed 
MACS samples is shown in Fig. [5] The histograms have been 
normalized to the number of clusters in the observed MACS 
sample (12). A two-sided Kolmogorov-Smirnov test returns a 
probability of ~ 12% that the observed and simulated data are 
drawn from the same distribution. The most remarkable differ- 
ence between the two samples is the excess of clusters with 
large lensing cross section among the MACS clusters. Three 
clusters, namely MACSJ0717.5+3745, MACSJ0025.4-1222, 
and MACSJ2129. 4-0741, have lensing cross sections larger 
than any cluster in the MareNostrum Universe. The medians 
of the two distributions are cr mM gyiACS = 3.07 x 10~ 3 /j~ 2 Mpc 2 
and cr medMNU = 2.7 x 10~ 3 /i~ 2 Mpc 2 for the observed and 
for the simulated MACS samples, respectively. The lensing 
cross section of cluster MACS J091 1.2+1746 is a factor of ~ 2 
smaller than the smallest cross section among the simulated 
clusters. As shown in Fig. [4] this cluster lays below the pre- 
dicted log <x - log %med scaling relation, i.e. the lensing cross 
section is small given the size of the Einstein ring. 




-3.0 -2.5 
Log(a) [h" 2 



Fig. 5. The distributions of the strong lensing cross sections. 
The blue filled histogram shows the results for the simu- 
lated MACS sample constructed with clusters taken from the 
MareNostrum Universe. The red shaded histogram shows the 
same distribution but for the observed MACS sample. The ver- 
tical dot-dashed lines indicate the medians of the two distribu- 
tions. 



The distributions can be used to estimate the differences 
between the expected number of giant arcs produced by the 
two cluster samples (when they lens sources at redshift z s = 2). 
The number of arcs given by 



(14) 



;=1 



where n s is the number density of sources in the background of 
the clusters, cr,- is the lensing cross section of the i - th cluster 
in the sample, and the sum is extended to all the clusters in the 
sample. Using Eq. [14] we estimate that about a factor ~ 2 less 
arcs are expected from the MareNostrum clusters compared 
to what expected from the MACS clusters. This is far from 



the or der-of-magnitude difference found by Bartelmann et al 



(1998) between clusters simulated in the framework of the 
ACDM cosmological model and observations, reducing sub- 
stantially the size of the arc-statistics problem. This is not sur- 
prising since our simulations inclu de the effects of mergers 
(Tom et al . 2004; Fede li et al which could n ot be prop- 

erly taken into account by iBartelmann et al.1 (119981) due to the 
limited number of clusters in their sample and to coarse time 
resolution in their simula tions. Analy z ing th e MareNostrum 
Universe, we showed in Fed eli et al. I (1201 Oh that un-relaxed 
clusters contribute to ~ 70% of the optical depth for clusters 
at z > 0.5. 

5.4. Distributions of Einstein ring sizes 

In Fig.|6]we compare the distributions of the Einstein ring sizes 
for the simulated and the observed MACS samples. As found 
for the lensing cross sections, the distribution derived for the 
observed cluster sample is also characterized by an excess of 
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Fig. 6. The distributions of the Einstein ring sizes. The blue 
filled histogram shows the results for the simulated MACS 
sample constructed with clusters taken from the MareNostrum 
Universe. The red shaded histogram shows the same distri- 
bution but for the observed MACS sample. The vertical dot- 
dashed lines indicate the medians of the two distributions. 



clusters towards the large values of #e- The same clusters hav- 
ing larger cross sections compared to any simulated cluster, 
also have Einstein radii exceeding the maximum value found in 
the simulations. The Kolmogorov-Smirnov test returns a prob- 
ability of ~ 30% that the two datasets are drawn from the same 
statistical distribution. The Einstein ring medians of the ob- 
served and of the s imulated samples differ by ~ 25%, while 
Zitrin et alJ (120 1 Oah report a difference of ~ 40% between 
the observed and the theoretical distributions of Einstein radii. 
However, their theoretical estimates are based on analytic mod- 
els of galaxy clusters, which are described by means of axially 
symm etric lenses with NFW density profiles. [Meneghetti et al. 
(2003) showed that such simple lens models under-estimate 
the lensing cross sections for giant arcs of numerically simu- 
lated clusters, mainly because they miss several important fea- 
tures like asymmetries and substructures, which enhance sig- 
nificantly the lensing efficiency of galaxy clusters (see also 
Meneghetti et al.l 120071) . Given the strong correlation existing 
between lensing cross section and Einstein radius, shown in 
Fig-El it is not surprising that our numerically simulated galaxy 
clusters have larger Einstein radii. Thus, as the arc statistics 
problem, also the Einstein ring problem results to be alleviated 
on the basis of our results, although three clusters still do not 
have a counterpart in our ACDM simulation. 



5.5. Expected concentration bias 

Meneghetti et al. (2010a) showed that, due to the cluster triax- 
iality and to the orientation bias that affects the strong lens- 
ing cluster population, we should expect to measure substan- 
tially higher concentrations than expected in clusters show- 
ing m any strong lensing features (see also Qgu ri et al. 2005, 



Fig. 7. The probability density function of the ratios between 
2D and 3D concentrations for clusters with the same properties 
of those in the MACS sample. The vertical dot-dashed line in- 
dicates the median of the distribution. The dashed line shows 
the cumulative distribution. 

timate the bias which is expected to affect a MACS-like sam- 
ple of clusters. For each of the clusters in the simulated sam- 
ple we measure the concentration inferred from both the pro- 
jected and the three-dimensional mass distributions, C2d and 
C3d- These are measured by the surface mass density and the 
three-dimensional mass density with NFW models. The con- 
centration bias is quantified by means of the ratio between the 
two-dimensional and the three-dimensional concentrations. 

The distribution of C2d/c3d derived from the simulated 
MACS sample is shown in Fig. [7] As expected, a median bias 
of the order of ~ 1 1 % is found. However, the distribution is 
skewed towards the high values. On the basis of our simula- 
tions we expect that the concentration may be biased by up to 
100% for some of the clusters. The dashed line shows the cu- 
mulative probability distribution. For about ~ 20% of the clus- 
ters in the simulated sample we measure a 2D concentration 
which is > 40% higher than the 3D concentration. 

In Fig. [8] we show how the concentration bias depends on 
the lensing cross section. While the bias is very small, of the 
order of few percent, for clusters with cr < 10~ 3 /i~ 2 Mpc 2 , it 
becomes increasingly higher for the clusters with larger cross 
sections. This is a very important result for interpreting the 
high conce ntrations recently measured in clusters like A1689 
or CL0024 (Broadh ursTetal]l2008l:lz"itrin et al.ll2009bh . 



2009). Indeed, the amplitude of this bias is a growing func- 
tion of the lensing cross section. It is worth (and easy) to es- 



6. Discussion and conclusions 

In this paper, we have used a large set of numerically simu- 
lated clusters, for which both the X-ray and the lensing proper- 
ties were previously investigated, and we have compared their 
Einstein ring sizes and strong lensing cross sections to those 
of clusters in an X-ray selected sample at redshift z > 0.5 
(MACS). We have also used the numerical simulations to 
estimate the expected bias in the lensing-inferred concentra- 
tions of the same MACS sample. Our aim was to investigate 
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Fig. 8. The concentration bias vs the strong lensing cross sec- 
tion. Each square represents a simulated cluster from a mock 
MACS catalog. The solid line is the median relation, while the 
error-bars indicate the inter-quartile ranges in each bin. 

whether some previously claimed discrepancies between the 
strong lensing properties of real and analytically modeled clus- 
ters in the framework of ACDM cosmology were confirmed by 
adopting more sophisticated and realistic cluster mass distribu- 
tions. These discrepancies between theory and observations re- 
gard 1) the large observed number of gravitational arcs behind 
the cores of galaxy clusters, 2) the large size of the Einstein 
rings reconstructed through lensing-based observations, and 3) 
the large mass concentrations determined by fitting the two- 
dimensional mass distribution of strong lensing clusters. Our 
main results can be summarized as follows: 

- the comparison between the numerical and the observa- 
tional cluster sample shows that some real clusters have 
too large lensing cross sections and Einstein rings com- 
pared to expectations in a ACDM cosmological model. 
However, the discrepancy between observations and sim- 
ulations is now significantly reduced compared to previous 
studies based on analytical estimates. Using N-body and 
hydrodynamical simulations allows to better capture sev- 
eral important properties of strong lensing clusters, such as 
their triaxiality, asymmetries, concentration scatter, dynam- 
ical activity, all of which were proven to boost the lensing 
efficiency; 

- using numerically simulated clusters we can predict the 
lensing concentration bias expected for the MACS sample. 
Such an estimate is based on the assumption that the numer- 
ical models are representative of the real cluster population. 
As shown by our results on the lensing cross sections and 
on the size of the Einstein rings, this may not be the case. 
Measuring a larger-than-expected concentration bias in real 
clusters may be a further indication that some substantial 
difference exists between the inner structure of observed 
and simulated clusters. 

We believe that these results are important because they 
suggest that our simplified models of the physical processes 



affecting the matter distribution in the cores of galaxy clus- 
ters must be improved. First of all this improvement concerns 
the description of gas-dynamical processes in clusters. On the 
other hand, our results may also indicate that some fundamen- 
tal assumption in the ACDM model is incorrect and leads to 
wrong predictions of the matter distribution on scales probed 
by galaxy clusters. In this case, alternative cosmologies, in- 
volving dynamical dark-energy, modified gravity, or primor- 
dial non-gaussianity may better explain the differences between 
the strong lensing properties of simulated and observed galaxy 
clusters. 

However, a few aspects of our work deserve some fur- 
ther discussion and deeper investigation. The simulations used 
in this study are performed in the framework of a WMAP- 
1 normalized cosmology. The value of <t& used in this study 
is larger than measured in the WMAP-7 data release. If the 
WMAP-7 normalization was adopted, the differences between 
simulations and ob servations would be significantly amplified. 
iFedeli et al.1 (|2008) show that, decreasing <x 8 from 0.9 to 0.8, 
the average lensing optical depth drops by at least a factor 
of three. The cluster efficiency for strong lensing would be 
strongly dimmed at z > 0.5, as a result of the delayed struc- 
ture formation caused by a lower normalization of the CDM 
power-spectrum. 

Our simulations suffer of two important limitations. The 
first is certainly given by the size of the cosmological box. Due 
to the relatively small box size, we may be under-predicting the 
number of very high-mass clusters in the sample. Moreover, we 
may be missing some systems that are particularly dynamically 
active. Consequently, our numerical sample may miss some of 
the most powerful lenses. We will investigate this issue by us- 
ing a sample of massive clusters extracted from a larger cosmo- 
logical box. The second limitation is due to the simple descrip- 
tion of the gas physics in the simulation, where no radiative and 
feedback processes are con sidered. On the basis of the recent 
results o f lMeadetalJd2010b . we are confident that this does not 
affect significantly the lensing properties of the clusters. It has 
been shown that the energy feedback from AGNs and super- 
novae counter-act and compensate the effects of cooling on the 
strong lensing cross sections. However, as we explained earlier, 
the X-ray scaling relations are strongly affected by such phys- 
ical processes. We introduced a correction to the X-ray fluxes 
based on the comparison between simulations and observations 
for accounting for this poor description of the gas physics. In 
absence of more sophisticated cluster models, this was the only 
viable approach, which is however based on several approxi- 
mations. 

As mentioned in Sect. [3] most of the multiple image sys- 
tems detected behind the MACS clusters do not have spectro- 
scopically confirmed redshifts. This lack of information makes 
these cluster mass reconstructions uncertain, although we be- 
lieve that our conclusions would be not affected dramatically 
(see the discussion in the Appendix lAV This situation will be 
soon improved thanks to the upcoming Hubble Multi-Cycle- 
Treasury-Program CLASljl, which will dedicate 524 new HST 
orbits to observe 25 galaxy clusters in 16 different bands, span- 



http : //http : //www . stsci . edu/~postman/CLASH/ 
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ning the near-UV to near-IR (P.I. Postman). These observa- 
tions, combined with already existing data for some of the tar- 
gets, are expected to deliver ~ 40 new multiple image systems 
per cluster with photometric redshift determined with an accu- 
racy Az ~ 0.02(1 Half of the clusters in the MACS sample 
used in this paper ar e also in the CLASH target list. As shown 
bv Meneghetti et alJ (|2010b), these observations will allow to 
constrain the Einstein ring sizes and the projected mass distri- 
butions in the inner regions of clusters with accuracies at the 
percent level. We look forward to improving our comparison 
between observations and simulations using this new dataset 
as well as improved cosmological simulations. 

While our pape r was in the process of being refereed, 
Horesh et al. ( 201 lb posted a paper were a comparison be- 
tween the lensed ac statisti cs by simulated hal os taken from 
the Millennium simulation ( Springel et al. 2005) and by a sam- 
ple of X-ray selected clusters is shown. Their results for z > 0.5 
agree with ours, being the simulated arc production efficiency 
lower by a factor of 3 than observed in the MACS cluster sam- 
ple. Instead, at lower redshift (z ~ 0.4), they find a very good 
agreement between the observed and the simulated arc statis- 
tics, in terms of the mean number of arcs per cluster, the dis- 
tribution of number of arcs per clusters, and the angular sepa- 
ration distribution. At even lower redshift, {z ~ 0.2) they again 
find an excess of arcs in the observations compared to simula- 
tions. Clearly, this emphasises that we need much larger sam- 
ples before we arrive at a firm conclusion. 
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Appendix A: How stable are the MACS lens 
models? 

As we pointed out earlier, the major source of uncertainty in the mass 
models of the MACS clusters used in this study is the lack of spec- 
troscopic redshifts available for most of the multiple image systems. 
Indeed, for several of these systems the redshifts are estimated photo- 
metrically. Other groups who performed a strong lensing analyses on 
some of the clusters used in this study did not agree on the identifi- 
cation of several strong lensing systems. In order to have an idea of 
how stable the models are with respe ct to the assumption s made, we 
compare here some of the models of Izitrin et al . (2010a) (Z models 
in the following) with alternative models kindly provided by Marceu 
Limousin. It is important to note that 1) these models are obtained with 
a completely different mass reconstruction code, the public software 
l enstooj^, which fol lows a completely different approach compared 
to Izitrin et al. Il l2010al) . and 2) in some cases the associations and the 
number of multiple image systems differ significantly. 

The lenstool models (L models hereafter) were pro- 
vided for four of the MACS clusters under investigation. 
These are the clusters MACSJ0717.5+3745 (Limo usin et al., 
in p rep.), MACSJ1 149.5+22 23 (see details in ISmithetalJ 
2009), MACS J0454. 1 -03 00 teerciano Alba et all I2007L |2010|) . 
and MACSJ1423. 8+2404 ( Limousin et alJl2010h . We use these mod- 
els fully consistently to the analysis made on the Z reconstructions, in 
order to derive measurements of the Einstein rings and of the lensing 
cross sections. 

The L model of MACSJ0717 is based on a set of 15 multiple im- 
age systems, the v ast majority of whic h agree with the systems iden- 
tified and used by Izitrin et al] ( feOlOal) . Recently Limousin and col- 
laborators gathered spectroscopic redshifts for two of these systems, 
which have been used to calibrate their model. The resulting recon- 
struction has and Einstein ring mc( j = 61" and a lensing cross section 
cr = 1.24 x 10 2 /T 2 Mpc~ 2 . These values are only ~ 15% smaller than 
those of the corresponding Z model. 

In MACSJ1149, the most spectacular lensed system is a spiral 
g alaxy whose total magnifica tion is estimated to be ~ 200. According 
to lZitrin & Broadhursfl ( Hooll and lZitrin etld1ll2010al) (Z model), this 
galaxy is lensed into five images, while lSmith et al. ( 2009) refer to the 
fifth image as part of the fourth image. This of course is a matter of 
interpretation. The calibration of the Z model is based on the assump- 
tion that the redshift of the above mentioned spiral galaxy is ~ 1.5, 



which was then confirmed spectroscopically by ISmith et al.l j2009h . 
Additionally, the Z lens model itself has allowed to identify six ad- 
ditional multiple-image candidate systems, that are not used in the 
construction of the L model nor are plausible in this model. The mass 
reconstructions derived from the constraints used by the two groups 
correspondingly differ. In particular, the Z model has a shallower pro- 
jected density profile than the L model. Despite this, the two models 
result to be almost identical in terms of both Einstein ring and lensing 
cross section, with differences of order of only a few percent. 

Using the L model of MACSJ1423 we find an Einstein ring and 
a cross section that are respectively ~ 30% and ~ 25% smaller than 
those derived from the Z model. Strangely, for this cluster spectro- 
scopic redshifts exist for two sets of multiple images which are used 
to build both the L and the Z models, thus such a mismatch between 
the models was quite unexpected. However, by visually inspecting the 
critical lines of the two reconstructions, we notice that the Z model 
predicts a northward extension of the tangent ial critical line, su rround- 
ing a bright elliptical galaxy (see Fig. 19 of Zitrin et al. 2010a). Such 
extension is not present in the L model and contributes to increasing 
both 8 E and cr. It may be an artifact or not, depending on the still un- 
known weight of this galaxy, as there are no lensing constraints in that 
region of the lens plane which probe this portion of the critical curve. 

MACSJ0454 shows the largest discrepancy between L and Z mod- 
els. In this case, the Einstein ring size and the lensing cross section de- 
rived from the L model are extremely larger than those derived from 
the Z model (a factor of two and a factor of five, respectively). We do 
not fully understand the origin of this mismatch. The Z model is built 
using the same lensing constraints as the L model, including also the 
spectroscopic redshift of one system at z s = 2.9. The Z model is also 
validated by additional two systems of multiply-lensed images, which 
were identified thanks to the model itself. Therefore, we are confident 
that the Z model used in the analysis is better constrained than the L 
model. 

To summarize, we find a good agreement between Z and L mod- 
els in three out of four of the MACS clusters used in this comparison, 
while for the fourth cluster, we think that the Z model is better con- 
strained, supporting our results. Therefore, we are confident that out 
results are robust. In the cr-6 E plane, even the L models do not depart 
significantly from the relation fitted to the simulations. This is shown 
in Fig. IA.ll which shows with sticks and diamonds the new locations 
of the 4 MACS clusters, when they are modeled using lenstool, and 
their shift with respect to the corresponding Z models. These are over- 
laid to Fig.|4]for comparison. 



3 see http : / /www . oamp . f r/cosmology/lenstool/ 




Fig.A.l. The shift of the some of the MACS clusters in the 
cr - 9 E , when alternative models obtained with lenstool are 
used to derive the Einstein rings and the lensing cross sections. 



